Generalised definitions of certain 
functions and their uses 

Debasis Biswas^, S. Biswas^'^ 
and 

Asoke P. Chattopadhyay^ 
^Department of Physics, ^Department of Chemistry 
University of Kalyani, 
Kalyani 741235, 
India 

Mbiswas-chak@rediffmail.com, ^asoke@klyuniv.ernet.in 
^sbiswas@klyuniv. ernet . in 

February 7, 2008 

Abstract 

Generalised definitions of exponential, trigonometric sine and cosine 
and hyperbolic sine and cosine functions are given. In the lowest order, 
these functions correspond to ordinary exponential, trigonometric sine 
etc. Some of the properties of the generalised functions are discussed. 
Importance of these functions and their possible applications are also 
considered. 
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1 Introduction 



Ordinary differential equations (ODEs) are ubiquitous in various branches 

of physical and biological sciences, and are essential to our understanding of 
topics as diverse as stability theory, analysis of electrical signals and networks, 
chaos and nonlinear dynamics, radiative processes, etc.[l, 2, 3, 4, 5, 6, 7, 8, 9]. 
In this paper we try to generalize the definition of sine, cosine and hyperbolic 
functions as solutions of differential equation of the form i/xx + f{x)y — 0. For 
f{x) oc x, x"^, x^, a;~^ the equations can be brought to Airy type. Coulomb 
type. Harmonic oscillator type and Euler type equations having wide appli- 
cation in physical science problems. We show that many well known ODE 
solutions of physical science can also be defined in terms of these generalized 
functions which we call now as g functions. 

The organization of the paper is as follows. In section 2 we formulate the 
definition of g functions related to generalised exponential, generalised cosine 
and sine and generalised cosine hyperbolic and sine hyperbolic functions. In 
section 3 we discuss some important ODEs of mathematical physics in terms 
of g functions. In section 4 we end up with a concluding discussion. 

2 Formulation of g function 

We begin with formal definitions of exponential, trigonometric and hyer- 
bolic functions as follows. The exponential function may be defined as solution 
of the differential equation 

Similarly, trigonometric sine and cosine functions may be defined as solutions 
to the following differential equation 

while hyperbolic sine and cosine functions are defined as solutions of 

<^ 
dx"^ 

Let us consider solutions to the following set of differential equations 

dy 
d 

Case I : Generalised exponential function 

Following Forsyth [5], we express gerenal solutions of (4) as 



where S denotes the operation of integration from to x i.e. Jq dx and P = x^. 
Then the solution of (4) becomes 

^(n+l) ^(2n+2) ™(3n+3) 

^~ ^ (n + 1) ^ (n + l)(2n + 2) ^ (n + l)(2n + 2)(3n + 3) ^ "' ^' 

The series on RHS of (8) is convergent for finite x and n. The only restriction 
on n is that n ^ —1. We can also solve equation (4) by separation of variables. 
The result is 

Iny = 



(n+1) 
so that 

y = e "+i (9) 

It is easily verified that the RHS of the two expressions, (8) and (9), are iden- 
tical for finite x and n. These can be used to define exponential function of 
order n as 

2,(3n+3) 

eJx) = 1 + + + + ... (10) 

' (n+1) (n + l)(2n + 2) (n + l)(2n + 2)(3n + 3) ^ ' 

Note that the ordinary exponential function is zeroth order of this generalised 
exponential function e„(x). Indeed, with n = 0, the above expression becomes 



eo(x)=e^ = l + - + - + - + 



However, the generalised exponential allows expressions for other values of n. 
Behaviour of en{x) for some values of n are shown in Figure 1. 
Case II : Generalised cosine and sine function 

Let us solve equation (5) with the same technique. We get two independent 
solutions as 

2/1 = 1 + ,52p + S'^PS^P + S^PS^PS^P + .... 
and 

y2^x + S^Px + S^PS^Px + S^PS^PS^Px + .... 

Here S = Jq dx as before and P — —x^. Then these solutions define gen- 
eralised cosine and sine functions for us as 



^(n+2) j.{2n+A) 
„n/ \ ^1 ± , ± , 

^^^^ (n + l)(n + 2) (n + l)(n + 2)(2n + 3)(2n + 4) ^ ^ 

2; (n+3) ^{2n+5) 

9si^) = 2; - (^^2)(n + 3) ^ (n + 2)(n + 3)(2n + 4)(2n + 5) ~ ^^'^^ 



these functions arc identical with ordinary trigonometric cosine and sine func- 
tions respectively. But expressions for other values of n are allowed. Behaviour 
of these functions for some values of n are shown in Figures 2 and 3 respectively. 

We now extend the above definitions for and gg{x) for negative val- 

ues of n as follows. g~"^{x) = xg^~^{^) leads to 

9c "(^) = 1 - _ 1)^^ _ 2) + _ i)(^ _ 2)(2n - 3)(2n - 4) ~ ^^^^ 

™(3-n) ™(5-2n) 

^ ' (n-2)(n-3) (n-2)(n-3)(2n-4)(2n-5) ^ ' 

Note that these expressions of g~'^ are also trivially obtained by substituting 
-n for n in corresponding expressions for g"'. It must be noted that g^^ arc not 
defined for all negative values of n, while g~"'{x) is not defined for n =1 and 
2, gs^ix) is not defined for n = 2 and 3. These functions are plotted for some 
negative n values in Figures 4 and 5. 
Case III: Generalised hyperbolic function 

For P — x", all the signs in the two series above becomes positive, and we 
get as solutions of equation (6), 



2; (n+2) 

9hc{x) = 1 + (^+i)(^ + 2) ^ (n + l)(n + 2)(2n + 3)(2n + 4) ^ ^^^-^ 

^(n+3) ^(2n+5) 

9hs{^) = ^ + ^ 2)(n + 3) ^ (n + 2)(n + 3)(2n + 4)(2n + 5) ^ ^^^^ 

where we have identified these solutions as generalised hyperbolic cosine and 
sine functions respectively. Note that these functions are also not defined for 
n < integers. It is trivial to show that ghd^) = cosh{x) and glsi^) = sinh{x). 
ghci^) is plotted for some values of n in Figure 6. The behaviour of g^^^lx), for 
the same n values, is very similar and is therefore not shown separately. 

We extend these definitions for generalised hyperbolic functions to nega- 
tive n values in the same manner as was done for generalised trigonometric 
functions. In other words 



9hc{x) = 1 + ^n-l){n-2) + (n - l)(n - 2)(2n - 3)(2n - 4) ^ ^^^^ 

^(3-n) ^(5-2n) 

9hs W = ^ + _ 2)(n -3)^ {n-2){n- 3)(2n - 4)(2n - 5) + - ^ 



define the generalised hyperbolic functions for negative values of n. As before. 



case of g^"' and g^^'' apply for these functions as well i.e. n 7^ 1,2 for gi,J^{x) 
and n 7^ 2,3 for g^^^^x). Behaviour of these functions for some negative n 
values are shown in Figures 7 and 8. 

It is clear that the above generalised definitions hold for positive and neg- 
ative real values of n (except those mentioned above). Solutions to the second 
order ODEs with —x^y on RHS are oscillatory. The oscillatory nature of these 
generalised functions is tested as follows. Leighton's oscillatory theorem[10] 
asserts that a differential equation 



{p{x)y') + q{x)y = 



is oscillatory in (0, 00) provided /q°° '^^^^ ~^ ^ /o°° q{x)dx — > 00. As both 
these conditions are satisfied in this case, the ODEs in (5) and (6) and hence 
their solutions are oscillatory. 

It can be proved that g^ix) and g^ix), as also g^d^) and g^si^) arc mutu- 
ally orthogonal in {—00, 00). This shall be dealt with in a later communication. 
We now propose that the functions e„(a;), g^ix), g]!^{x), g^d^) and g^si^) be 
considered as generalised exponential, generalised cosine and sine, and gener- 
alised hyperbolic cosine and sine functions respectively. Prom plots of these 
functions in Figures 1 to 5, identity of these functions with ordinary exponen- 
tial, trigonometric cosine and sine, and hyperbolic cosine and sine functions 
for n = 0, and departure from the latter functions for higher values of n are 
evident. 

3 Applications 

(I) We next show the utility of these functions. Consider any second order 
ODE of the form 

d'^u du dR ^ , , 

— + 2R—^{R^ + — + x'')u^Q 19 
dx'^ dx dx 

Substituting u = ye~-^ , R=R(x), the above ODE becomes identical with 
equation (5), whose solutions (except for negative integral n) can be written in 
terms of g^{x) and g^{x). Therefore (except for some n < o integers), solutions 
of (19) can be expressed in terms of the latter functions. 

(II) These oscillatory functions can also serve as solutions to the following 
type of ODEs 

where f{x) = Ax"^ + Bx''^+ .... where A, B, .. are constants and m, n, ... may 
be relatively coprime. 



Mathieu's etc. We discuss briefly on these connections, details will be placed 
in future works, with applications of g^^{x) and gg{x) wherever possible. Here, 
we shall show how these functions can be used as solutions of Bessel's equation. 
Bessel's ODE is written as 

X -j-^ + X— + {x -u )y^0 
ax^ ax 

Putting X — 2vr^ and y — ur~^, the above ODE becomes 

(fu _j_ 

+ r—^u = 

the solution of which can be written as 

u{r) = Ag~'{r) + Bgr^{r) 



Hence, the general solution of Bessel's ODE can be expressed as 
y(x) = {^-^nAgr'i^r + Bgr'i^r] 



(IV) Let us discuss the solution of Ricatti equation in terms of g function. 
The Ricatti equation is 

+ y2 + x"* = (20) 



dx 

With — = I/, the above equation can be cast into the form 



+ x"* = (21) 
u 

which evidently has two solutions in terms of g^ and g^. 
(V) The confluent hypergeometric equation is 

If we substitute x = —t'^/4:, y = t^^z with /c = n + 1, the confluent hypergeo- 
metric equation is reduced to Bessel differential equation 

dj^ z dz 

and can be solved in terms of g^ and g^ as shown earlier. 

Moreover Whittaker equation and Weber-Hermite equation can also be 
solved in terms of g^ and g^. The above examples merely shows the equivalence 



We now present an example where the present method is shown to be su- 
perior to existing ones. In the case of an ODE such as 

d'^y -4 n 



no regular solution can be obtained by the usual, Probenius method. Only 
normal integrals exist in such cases [5]. But the present method immediately 
provides regular solutions in a straightforward manner. This is also appar- 
ent in cases where logarithmic terms appear in the Frobcnius method which 
are troublesome to deal with. Such a possibility does not arise in the present 
method. Also, for comphcated nonlinear ODEs for example, often one solution 
appears in published tables of integrals; the present method yields all possible 
solutions. We just cite one example. 

Consider the equation 

yy"' + yy" + x^yy'^Q (24) 
With substitution w — yy', the above equation reduces to the form 

w" + x'^w = (25) 
Obviously the solutions are Wi = g^{x) and W2 — 9^{x). 

4 Concluding discussion 

Applicability of the generalised exponential and the generalised hyperbolic 
functions defined above can also be established in a similar manner. It be- 
comes apparent that such functions can be applied to a fairly broad category 
of physical problems. While the present work may be considered as a brief 
introduction to these functions, their detailed nature and properties will be 
taken up in our future works. We hope to establish the usefulness of these 
functions, and show how solutions of various the problems of mathematical 
physics can be expressed in much more compact and uniform manner with 
their help. 



We wish to emphasise that the solutions to a class of differential equations 
given in the present work is actually a subset of a larger class of functions. 
Here, we have only considered first and second order ODEs of a particular 
type. We wish to state that these are subsets of the following ODEs viz. 

where n and m may both vary. We wish to express solutions to such ODEs as 
9nmi^)^ where k = 0, 1, ... n. For example, e„(a;) is now written as gifi{x), 



adequate for denoting solutions to lower order ODEs as shown at the beginning 
of this work, viz. equations (4) to (6), the need for a more general notation 
becomes immediately apparent as higher order ODEs are considered. For ex- 
ample, for 

the solutions can be written as g^:i^{x), g^m{x) and g'^J^{x). Such functions 
will be considered in our future publications. 
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Figure 1: e„(x) vs x for n = (red), 2 (green) and 10 (blue). 




Figure 2: g^ix) vs x for n = (red), 2 (green) and 10 (blue). Note increased 
oscillation as n is increased. The abrupt termination for the blue curve is due 
to numerical instability. 
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Figure 3: (/"(x) vs x for n = (red), 2 (green) and 10 (blue). Note increased 
oscillation as n is increased. The abrupt termination for the blue curve is due 
to numerical instability. 
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Figure 4: g^{x) vs x for n = -4 (red) and -10 (green). Note oscillatory be- 
haviour for small x and asymptotic behaviour for large x as expected from 
equation (13). Oscillations increases as magnitude of n is increased. The 
abrupt termination, more apparent for n = -10 is due to numerical instability. 
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Figure 5: g^{x) vs x for n = -4 (red) and -10 (green). Note that for small x, 
gs{x) ~ X as expected from equation (14). For larger x, the function oscillates, 
more for n = -10. The abrupt termination is due to numerical instability. 
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Figure 6: g^^ix) vs x for n = (red), 2 (green) and 10 (blue). 
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Figure 7: g^d^) vs x for n = -4 (red) and -10 (green). Note the instability 
near x = and smooth approach towards unity (first term in the series) for 
large x, as expected. 
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Figure 8: g^si^) vs x for n = -4 (red) and -10 (green). Note the instability 
near x = 0. For large x, g][g{x) ~ x as expected (see also Figure 5). 



